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Abstract 

Ab initio calculations of relativistic electron energy loss spectra (REELS) are carried out using 
a generalization of the real-space Green's function code FEFF8 which is applicable to general 
aperiodic materials. Our approach incorporates relativistic effects in terms of the cross-section 
tensor within the dipole selection rule. In particular the approach explains relativistic corrections 
to the magic angle in polarized EELS experiments. Our generalization includes instrumental effects, 
such as the integral of the cross section over the impulse transfer dependence, to account for a finite 
detector aperture and electron beam width. The approach is illustrated with an application to the 
graphite C K edge. 
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I. INTRODUCTION 



Electron energy loss spectroscopy (EELS) measures the energy loss of a fine beam of high 
energy electrons (~ 100 keV) propagated through a sample in an electron microscope.- The 
energy loss spectrum is defined as the fraction of electrons which has lost a given amount 
of energy by interacting inelastically with the sample. From the EELS spectrum, one can 
obtain structural, chemical and electronic information similar to that in x-ray absorption 
spectroscopy (XAS), encompassing both extended x-ray absorption fine structure (EXAFS) 
and x-ray absorption near edge structure (XANES). We focus here on the ELNES (energy 
loss near edge structure) edges in the spectrum, corresponding to inelastic losses through the 
excitation of an electron from a deep core level into unoccupied states. Because EELS is an 
absorption technique, and the initial core level state is sharp in energy, such core loss signals 
reflect the electronic structure of unoccupied electron states. In particular, apart from a 
smooth background cross-section factor, the shape of the ionisation edge is roughly an image 
of the unoccupied angular momentum projected density of electron states (LDOS). Due to 
selection rules, the observed spectra correspond to a decomposition of the DOS according 
into various angular momentum / channels. Also the possibility of tilting the specimen with 
respect to the electron beam at fixed scattering angle allows one to investigate the anisotropy 
in the local unoccupied DOS because the scattering vector q (i.e., the momentum transfer) 
appears in the transition matrix element. Here q = k — k', where k is the wave vector of 
the incident fast electron and k' the wave vector of the scattered fast electron. This is the 
analog of the linear dichroism well known in x-ray absorption spectroscopy (XAS). 

What distinguishes EELS from XAS and similar techniques is that one can now obtain 
very local atomic scale information by focussing a very small probe of width ^0.1 nm on a 
sample in a transmission electron microscope (TEM). In the age of nanotechnology, this is 
relevant for studies of nanoscale materials. Modern instruments with field emitters also allow 
detection of ELNES with an energy resolution of 0.6 to 0.7 eV on a sub-nanometer scale, 
and monochromated transmission electron microscopes (TEM) reach 0.1 eV resolution.- 

Ab initio calculations of EELS have often been used to support the interpretation of 
experimental data. Several approaches for these calculations have been developed. For 
example, one approach for periodic structures is based on density functional theory (DFT) 
and the LAPW band-structure code WIEN2K^'^ and super-cell techniques. An alternative 
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approach for EELS makes use of the real space Green's function based FEFF8^, which is 
apphcable to periodic and non-periodic structures ahke. Although this ab initio code has 
been used extensively in the field of x-ray spectroscopy, it has been frequently applied to 
EELS calculations as well. A recent review of the application of FEFF8 to EELS calculations 
can be found in Ref. [6;. This approach is based on the approximate equivalence between 
dipole-selected EELS (i.e., the long wavelength limit) and XAS that has long been taken 
for granted. In this paper, however, we examine the quantitative differences between XAS 
and EELS. In particular we discuss the differences due to relativistic effects as well as 
instrumental effects (e.g., characteristics of the electron microscope). These are needed 
to obtain quantitative agreement with modern EELS experiments using relativistic beam 
energies. 

Recently, it has been recognized*^ that a relativistic interaction Hamiltonian is essential 
for accurate calculations of the scattering cross section, i.e., to obtain quantitative results for 
anisotropic materials at relativistic beam energies. The effect is a relativistic compression 
of the interaction field, which is therefore anisotropic in the dipole limit. Thus, treating this 
relativistic effect requires a generalization of the FEFF8 code to account for the momentum- 
transfer dependence of relativistic EELS experiments. The main purpose of this work is 
to develop an approach for ab initio EELS calculations that builds in the relativistical 

u 

formalism of Ref. |7| within the framework of Green's function multiple-scattering theory. 
In particular this effort builds on the self-consistent electronic-structure/spectroscopy code 
FEFF8,- but introduces a number of extensions relevant to modern EELS experiments. This 
development is broadly applicable and provides a new, general tool for EELS calculations 
which is complementary to band-structure techniques and is applicable over a broad spectral 
range. 

One important application is the so-called magic angle, which has been at the heart of 
the discussion in the literature on relativistic EELS.- This angle is defined as the value of 
the detector aperture of an electron microscope for which the measured EELS spectrum is 
independent of the relative orientation of sample and electron beam. This quantity is of 
direct practical importance for polarized EELS experiments, in which the anisotropy of the 
signal can be an unwelcome complication. Notably the magic angle is material independent, 
and therefore provides a direct test of the validity of scattering theory that our calculations 
are based on. We show in this paper that relativistic calculations based on the new FEFF8 
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code significantly improve on nonrelativistic calculations. 

The theory relevant to the relativistic extension of FEFF8 is briefly described in Sec. [Tll 
while Sec. HII] describes computational considerations. Finally applications are given in Sec. 



II. THEORY 



The EELS signal can be described by the double differential scattering cross section 
(DDCS) 

(q, E) 

dndE ' ^ ^ 

which is the probability of detecting a scattered electron which has lost energy E and 
transferred-momentum q into the solid-angle dQ. Formally the DDCS can be expressed 
in terms of the bare Thomson cross-section and the relativistic dynamic structure factor 
S'(q, E). Since the Thomson cross-section is sharply peaked at small q, it is common practise 
and generally a good approximation to consider only the so-called dipole transitions (i.e., 
small q limit) where the orbital momentum quantum number / of the atomic electron changes 
by ±1 in transitions. Recently it has been shown'' that within the dipole approximation the 
relativistic DDSCS for EELS is given by 



d^a{q,E) _ fda_' 
dEdVL ~ \dVL 



S(q',E) (2) 



Th 



(da_\ ^ Aaf^^ k' 

\dnj^^ [q^-{E/hcffk ^' 

S(q',i?) = Y.m^-Am{Ef-E,-E). (4) 

Here the momentum transfer in the dipole transition element is relativistically contracted, 
i.e., 

q' = q - I3'^qz^2. (5) 

where (3 = v/c and v is the beam velocity. This equation is very similar to the description 
of XAS in the dipole limit, with the impulse transfer q playing the role of the polarization 
vector e in x-ray scattering matrix elements. However, for relativistic EELS there is an extra 
q-dependent contribution along the direction of propagation e^. 



In general the DDCS can always be separated into a probe-dependent part containing the 
q-dependence, and a sample-dependent part which is independent of q. Since the theory is 
bi-linear in q the sample-dependent term transforms as a tensor, i.e., 

S(q',E)= E q[q',a,,{E) (6) 

a,,{E) = (i\ X, I/) (/I X, |z) 5{E - E, + Ej) . (7) 

The cross-section tensor aij (CST) can therefore describe all possible transitions of the 
sample. However, experimental conditions determine which impulse transfers occur and 
therefore the weight of each component of the cross-section tensor that contributes to the 
total cross section. This can be illustrated clearly by considering the sample to beam ori- 
entation of an EELS experiment. Rotation of the sample is equivalent to a rotation of q, 
thus changing the weights of the aij components in Eq. ([6]). The relativistic character of 
the formalism is also obvious: the field of the beam electron contracts in its propagation 
direction, resulting in the evaluation of M using a contracted impulse transfer vector as in 
Eq. ([5]), denoted by a prime. 

Formally the CST is a symmetric tensor with at most six independent components. 
As such, it can always be diagonalized. However, only in symmetric materials, where its 
principal axes are given by the physical symmetry of the crystal itself, a priori knowledge of 
the diagonal representation is available, and one set of coordinates diagonalizes the tensor 
for all energies. In the general case of a low symmetry sample, or in a situation where 
a non-symmetric coordinate system is desirable, the cross terms in Eq. are important 
contributions to the cross section which cannot be neglected. We give an example of this in 
Sec. m 

III. PRACTICAL CONSIDERATIONS 
A. Description of the experiment 

An ab initio EELS calculation requires not only a description of the sample as input. It 
also requires knowledge of the conditions in which the experiment was done. We briefly dis- 
cuss some experimental parameters specific to EELS. An EELS experiment usually involves 
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an integration over the DDCS defined in Eq. ([2]). Typically, the probe has a certain angular 
width, characterized by the convergence semi-angle a, allowing a whole set of ki incoming 
plane waves. Similarly, the detector integrates the signal over a certain range of outgoing 
beam directions, characterized by the collection semi-angle (3. Both are usually of the order 
of mrad. Assuming that the incoming beam is monochromatic, the measured signal is then 
given by 

[q,E)d'q (8) 



dE J dVLdE 

a;l3 

. To cast the orientation dependence of the EELS spectrum in a more explicit form, we can 
rewrite Eq. ([8]) using the CST : 



- 1^^ ^^^(E)J^ yr^^^j^.^d q (9) 

As the CST depends only on the sample, only integrals of functions of q need to be 
calculated. These integrals are approximated by a sum over a finite set of impulse transfer 
vectors q. This is considered in more detail in the Appendix. 

Other experimental parameters included in our calculations are the electron beam energy, 
the sample to beam orientation, and the position of the EELS detector in the scattering 
plane. FEFF8 calculations always include core hole broadening. Additional broadening can 
be applied. 



B. Computational details 

Our calculations are incorporated in a generalization of the ah initio real space multiple 
scattering code FEFF8, which we call FEFF8.5. We made modifications to FEFF8 to obtain 
the full-cross section tensor cTjj from the program, and added a new module that calculates 
the net cross section of Eq. ^ for a given q from the cross section tensor, and carries out the 
integrals over q as described in Sec. IIII A[ We have taken care to retain all other features 
of the code, including the use of advanced cards such as TDLDA in the calculations^ which 
account for corrections to the independent electron approximation ; the use of Debye- Waller 
factors to approximately account for temperature effects ; etc. 

The calculation of the cross section tensor for EELS is analogous to the case of XAS 
calculations,^ For the near edge region or energy loss near edge structure (ELNES) - 
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the electron equivalent of XANES), the full multiple scattering technique (FMS) is used, 
in which all scattering paths within a sphere of limited radius are summed implicitly by 
matrix inversion. For the extended region or extended energy loss fine structure (EXELFS) 
- the electron equivalent of EXAFS, the path expansion approach is taken, in which the 
scattering from a selected number of paths of limited length is summed explicitly. This is 
done for each of the six independent components of the sigma tensor. Combining the ELNES 
and EXELFS calculations, one can calculate spectra over hundreds of eV, far beyond the 
limitations of most band-structure codes. 

Note that the calculation of the probe and the sample are treated separately so that 
it is sufficient to calculate the properties of the sample once for the simulation of many 
experimental situations. We remark that our EELS code can also calculate mixed dynamic 
form factors (MDFF), which are off diagonal in q, but we postpone the discussion of these 
to a future publication. 

IV. APPLICATIONS 

A. The C K edge of graphite 

The strongly anisotropic nature of graphite makes its EELS spectra highly susceptible 
to orientation effects and therefore to relativistic contributions to the cross section that we 
discussed in Sec. HIl We demonstrate our method on the C K edge, which has an energy 
threshold of 285 eV. 

1. Cross-section tensor 

In this subsection, we discuss the components of the CST, and show that calculation of 
its diagonal components is not generally sufficient to calculate the EELS spectrum. To see 
this, we show the different components of the cross section tensor calculated in two different 
coordinate systems. System 1 is symmetrical : its z-axis is perpendicular to the graphene 
sheets of the sample, x and y are in-plane. System 2 is non-symmetrical : it is obtained 
from system 1 by a rotation of 35 around the x-axis of system 1. We work in a Carthesian 
representation and refer to the components i,j of a as x,y,z. 
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graphite symmetric frame - XANES sigma tt 
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FIG. 1: Components of the CST of graphite in symmetric coordinates. 
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FIG. 2: Components of the CST of graphite in nonsymmetric coordinates, tilted 35 around the 
symmetric x-axis. 



In Fig. [H we see that in symmetric coordinates the spectrum contains the so-called 
TT-transitions, while a^x and ayy are identical and contain the a-transitions. All off-diagonal 
components are zero, as can be explained by symmetry, i.e., equivalence of x and -x, y and 
-y, and z and -z. 

Fig. [2] shows that in the rotated system, a^x is equal to that in the symmetric frame, 
but ayy and azz have mixed and are of mixed vr, cr - character. Additionally, the decrease 
in symmetry allows y,z cross terms to exist. The x, -x symmetry has been preserved, 
suppressing xz, zx, xy and yx components. A more general rotation of the coordinates 
would make all off-diagonal elements nonzero. 

Finally, Fig. [3] shows the resulting ELNES spectrum. In system 1, calculation of the direct 
components of a is sufficient. To calculate the same spectrum in system 2, however, the off 
diagonal components {yz and zy in this example) make a very important contribution. 



8 



graphite C K edge 











symmetric frarr 


s — 1«— - 












rotated frame 


cross tetui 


S B 
























1 
































\ 














/ 


















3BBEQEini3 


l<'« II II II'M'II' 








I 















270 280 290 300 310 320 330 340 350 

energy loss in eV 



FIG. 3: C K edge of graphite calculated in coordinate systems 1 and 2 (see text). The beam is 
perpendicular to the graphene sheets, the beam energy is 300 keV, q = 10 mrad, /? = mrad. 

2. The Magic Angle 

The magic angle is defined as that value the collection angle jSm of for which the EELS 
spectrum is independent of sample to beam orientation, given a certain convergence angle 
a. In the dipole approximation, one can prove easily that such an angle exists^ at which the 
integrals in Eq. (Q lose their orientation dependence. The magic angle depends only on beam 
energy and energy loss (but it is approximately constant over the near edge region) and not on 
any material property. The magic angle has played a key role in recent developments of EELS 
theory^ - arguably exactly for that reason. It is of practical importance to experimentators 
wishing to eliminate the complications of orientation dependence from their investigations, 
and turns out to be sensitive to the details of scattering theory. It is often expressed in units 
of the "characteristic scattering angle" 6e, which is the width of the Lorentzian function 
that approximately gives the DFF as a function of scattering angle 

where Eq is the beam energy and rrie is the electron rest mass. In Fig. HI the magic angle 
is experimentally found to be /3m = 0.68 mrad for given experimental conditions.— This is 
close to 9e = 0.59 mrad. 

We now turn to theoretical calculations using FEFF8. We calculate spectra at different 
sample to beam orientations, which we characterize by a single tilt angle between the electron 
beam and the crystal c-axis. This tilt corresponds to a rotation of q'j in Eq. [9l We could 
investigate the rotation invariance of differential cross section, but it is more convenient to 
choose a more sensitive function of the spectrum, and study it as a function of collection 
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FIG. 4: C K edge of graphite measured^- at 300 keV beam energy for 2 orientations, and different 
collection angles. The two orientations overlap between the camera lengths of 80 cm and 100 cm, 
corresponding to a magic value of the collection angle of about 0.68 mrad. The convergence angle 
is close to mrad. 

angle at fixed energy loss. At the magic angle, the partial i = j cross sections of Eq. [9] 
are individually rotation invariant,^ and therefore we may equivalently study the ^ ratio of 
the spectrum. This function is important as it is related to the sp"^/ sp'^-ratio which is often 
used to characterize carbon samples.— In symmetric coordinates, it is given by the vr term 
or i = j = z term of Eq. divided by the a term or i = j = x plus i = j = y term. 



1-1 



(11) 



We calculate this quantity at a fixed energy loss oi E = 29AeV as a function of collection 
angle (3, shown in Fig. O for a nonrelativistic calculation, and in Fig. E] for a relativistic 
calculation. Three different sample to beam orientations are shown in each Figure. At the 
magic angle, the spectrum and its ^ ratio are independent of orientation. 

The nonrelativistic simulation in Fig. [5] gives = 4=0 e, as has been reported in the 
literature for many yeara^, but is inconsistent with experiment. The relativistic calculation 
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FIG. 5: ^ ratio (see text) of the graphite C K edge at 10 eV above threshold for a 300 keV beam 
and three sample to beam orientations. The magic angle is at the intersection of the three curves. 
Nonrelativistic calculation. The convergence angle is mrad. 
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FIG. 6: ^ ratio (see text) of the graphite C K edge at 10 eV above threshold for a 300 keV beam 
and three sample to beam orientations. The magic angle is at the intersection of the three curves. 
Relativistic calculation. The convergence angle is mrad. 

in Fig. E] yields a magic angle f3rn = 0.65 mrad, in much better agreement with experimental 
measurements giving about 0.68 mrad in Fig. HI The same information is contained in figs. 
[7] and [HI where relativistic calculations of the spectrum are shown at both values of the 
collection angle. The nonrelativistic calculation would yield identical spectra at /9 = 2.4 
mrad, in disagreement with experiment. 

Our present results agree very well with calculations reported in^^, which were calculated 
using the DFT code WIEN2K. 
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FIG. 7: C K edge of graphite for 3 orientations, 300 keV beam energy, q = mrad, /? = 2.4 mrad). 
Relativistic calculation. 
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FIG. 8: C K edge of graphite for 3 orientations, 300 keV beam energy, q = mrad, /? = 0.6 mrad). 
Relativistic calculation. 

V. CONCLUSIONS 



We have presented relativistic calculations of electron energy loss spectroscopy using the 
real space Green's function code FEFF. The calculations correctly accounts for q-dependence 
and microscope settings such as collection and convergence angle. We have demonstrated 
our method on the C K edges of graphite, where we calculate the correct magic angle. 
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APPENDIX A: INTEGRATING THE CROSS SECTION OVER BEAM AND 
DETECTOR 

In this section, we consider the calculation of the differential cross section (Eq. ([8]) in 
more detail. For brevity, we will write d{E) or d{k') for the differential cross section, and 
d'^{E,Q) or (i^(k, k') for the double differential cross section, where k and k' are the wave 
vectors of the incoming and outgoing beam electron, and q = k — k' is their difference. The 
detector opening is a circle of radius (3 and the electron beam that hits the sample has a 
profile fa- The differential cross section d{E) = d{k') is then given by 

d{E) = J dQ' J d^kd\k,k') fa{k) (Al) 
= Jd^d'ici) J^fa{k' + q)dn' (A2) 
= Jd'qd\q)g'^^{k',ci). (A3) 

We have used the fact that c?^ is invariant to small rotations of the order of typical scattering 
angles in EELS (of order mrad), and hence depends only on q. A uniform, monochromatic 
(of fixed energy /cq), circular beam is described by 

Uk) = ^0^(9) 5ik - ko) (A4) 

With this type of beam, the weight g in Eq. (lAip is an integral over a constant function, 
i.e., a surface. If the detector aperture is a circle of radius /?, the weight can be interpreted 
as an overlap of two circles of radius a and (3 whose centers are separated by the vector qj_. 

g^^'ik', q) = A n {Qa - q±) (A5) 

which can readily be evaluated using basic algebra. If collection and convergence angle are 
interchanged, the shape of the spectrum is conserved, but it is multiplied by (a//9)^. More 
complex beam profiles (or detector aperture profiles) could destroy this pseudo-equivalence. 



Electronic address: kevin . j orissen@ua. ac . be 



13 



^ R. F. Egerton, Electron Energy-Loss Spectroscopy in the Electron Microscope (Plenum Press, 

233 Spring Street, New York, N.Y. 10013, 1996). 
^ S. Lazar, G. Botton, M.-Y. Wu, F. Tichelaar, and H. Zandbergen, Ultramicroscopy 96, 535 

(2003). 

3 C. Hebert, Micron 38, 12 (2007). 

^ P. Blaha, K. Schwarz, G. Madsen, D. Kvasicka, and J. Luitz, WIEN2k, An Augmented Plane 
Wave + Local Orbital Program for Calculating Crystal Properties (Karlheinz Schwarz, Techn. 
Universitat Wien, Austria, 2001), ISBN 3-9501031-1-2. 

5 A. L. Ankudinov, B. Ravel, J. J. Rehr, and S. D. Conradson, Phys. Rev. B 58, 7565 (1998). 

^ M. Moreno, K. Jorissen, and J. Rehr, Micron 38, 1 (2007). 

^ P. Schattschneider, C. Hebert, H. Franco, and B. Jouffrey, Phys. Rev. B 72, 045142 (2005). 

^ A. L. Ankudinov, Y. Takimoto, and J. J. Rehr, Phys. Rev. B 71, 165110 (2005). 

^ C. Hebert, P. Schattschneider, H. Franco, and B. Jouffrey, Ultramicroscopy 106, 1139 (2006). 

^ J. Verbeeck, University of Antwerp (unpublished). 

1 J. T. Titantah and D. Lamoen, Phys. Rev. B 70, 075115 (2004). 

^ K. Jorissen, J. Luitz, and C. Hebert, Ultramicroscopy (submitted). 

^ See the FEFF documentation |http://leonardo.phys.washington.edu/feff/ 



14 



